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The Identification Problem 



Measured Input u(t) Measured Output z(t) 
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Identification of Stochastic or Noise Model 



Measured Input u(t) Measured Output z(t) 


- Relationship between w\(t) and z(t), given only system output, z(t) 

- u(t ) assumed zero or constant 

- Commonly known as time-series analysis 

- Economic analysis, geophysical or astronomical phenomena, biological data (e.g. EKG, 
EEG), etc. 
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Identification of the Deterministic Model, Q 



Measured Input u(t) Measured Output z(t) 


- Relationship between u(t) and y(t) - assumes w\(t) = 0 

- Input/output corrupted by noise, e u (t) & e z (t) - commonly assumes e u (t) = 0 

- Pursued when objective is to gain insight into the functioning of a system 

- Automotive industry, chemical plants, pulp & paper, biomedical modelling (e.g. respira- 
tory modelling, drug delivery, disease modelling, etc.) 
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Identification of Stochastic & Deterministic Models 



Measured Input u(t) Measured Output z(t) 


- Both input/output signals available for identification 

- Used when accurate predictions are desired 

- Design of model-based control systems for aircraft, spacecraft or robotics. 
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Representation of the System Q 


• Linear or Non-linear 

- Systems in nature are inherently non-linear - especially in biology 

- Linear modelling about an operating point is difficult with biological systems due to 
operating point variability 

- Could lead to misinterpretation of physiology, e.g. disease quantification 

- Biological processes should be modelled in their natural state - non-linearly 

• Nonparametric or Parametric 
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Nonparametric Methods 

.7 — — ^ 7 -- i: ifc, ~ »« E 


• Advantage 

- Provide convenient, robust means of characterising the dynamics of linear systems 
without requiring a priori assumptions regarding the system structure 


• Disadvantage 

- Nonparametric estimates of dynamics are difficult to relate to the structure and param- 
eters of the underlying physiological system 
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Parametric Methods 


• Disadvantage 

- Generally require a priori assumptions about the system order 


• Advantage 

- Provides a concise description of the system dynamics 

- Yield results that may be related directly to the system structure 
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Model Form 


■JVA 


• Linear statistical model 

z (n) = ^2e 3 f(ip 3 (n)) + e(n) 

.7=1 

- 2 : observed system output 

- Qj : unknown system parameter 

- /: non-linear mapping 

- iff regressor 

- e: independent Gaussian variable, zero-mean, <r 2 homoskedastic (constant, finite vari- 
ance) 

• Let ip be described as: 

T 

<p{n) = [1, u(n), • • • ,u(n- n u ), z(n - 1), • • • ,z(n- n z ), e(n - 1), • • • , e(n - n e )] 

- Special case / polynomial: u 2 (n — 3), u(n)u(n— 1), z(n— l)z(n — 2), u 2 (n— l)z(n — 2) 

- General case /: wide variety of non-linear functions such as a sigmoid 

- NARMAX 

• Linear-in-the-parameters 

- Linear or pseudolinear regression techniques 
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Example of a NARMAX Model 


• Non-linear feedback model 


Flap Angle 





&1S 3 +&2<s 2 +&3S+64 

s 4 +ais 3 +a2S 2 +a3S+a4 


dX 2 {s) 




• NARMAX representation 


Noise-Free Wing Pitch Response 



Pitch Measurement Error (Noise) 



0 10 20 30 

Time (s) 


z(n) — 6iz(n — 1 ) + 62z(n — 2 ) + 6^z(n — 3 ) + 6^ z(n — 4 ) + 6^z 2 (n — 1 ) + 9 §z 2 (n — 2 ) 
+ 0 jz 2 {n — 3 ) + 6%z 2 (n — 4 ) + 6gz 2 (n — 5 ) + 6iou(n) + Q\\u(n — 1 ) + 9\2u(n — 2 ) 
+ 0 isu(n — 3 ) + 6uu{n — 4 )— 0 ie(n — 1 ) — 62 e(n — 2 ) — 6^e{n — 3 ) — 9 ^ e(n — 4 ) 

— 26^z(n — 1 )e(n — 1 ) — 26qz(ti — 2 )e(n — 2) — 26jz{n — 3 )e(n — 3 ) 

— 26%z(n — 4 )e(n — 4 ) — 26gz(n — 5 )e(n — 5 ) + 6$e 2 (n — 1 ) + 0Qe 2 (n — 2 ) 

+ 0^e 2 (n — 3 ) + 0%e 2 (n — 4 ) + 0ge 2 (n — 5 ) + e(n) 
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Full Identification Procedure 




• Model order selection 

- Determine number of input, output and error lags and non-linearity order/basis function 

• Parameter estimation 

- Determine values of unknown parameters 

• Structure detection 

- Select parameters to include in model 

• Model validation 

- Assess whether identified nominal model can reproduce data from a plant 
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Model Order 


z{n) = f[l,u(n), 


,u(n - n u ),z(n - 1), 
e(n - I),-' 


• , z(n — n z ), 
, e(n — n e )\ + e(n 


• Model order represented as 

O 


Ti u Tlz Tie l 
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Parameter Estimation 

Cfc— ^ - • I ■ ■ — ~ i: ifa E 


• Need an estimate of 0 using standard £ 2 minimisation 


1 

min - 
e 2 




2 

2 


• NARMAX models provide concise system representation 

- Noise on output enters model as product terms with system input and output making 
parameter estimation challenging 


Ordinary least-squares yields biased estimate: does not account for noise 


6* = (^ T ^) _1 ^ T Z where = [<!>,,, 


• Solution extended least-squares (ELS) 

6 = (’$' T ^) 1 ^ T Z where ^ 
- Bias addressed by modelling lagged errors 


^ zu & > & 


zue 
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Structure Detection 


• Selection of subset of candidate terms that best describe observed output 

• Maximum number of candidate terms 

i 

p = ^2pk + 1 

k = 1 

Pk-i(n z + n u + n e + k) 

Pk = 7 ■ Po = 1 

- Example: model of order: O = [4 4 4 2] 

- p = 105 candidate terms 

- The curse of dimensionality! 

• Leads to computationally intractable combinatorial optimisation problems 
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Example of Candidate Terms and Structure Detection 

’■ i * irr. I ■ ■ ■ ■ - -- M - ! tit -,jj B 


• Model example: 


z(ri) = Oiu{n — 1 ) + 9 l u 2 (n — 1 ) + 9 l z{n — 1 ) + 9 % e[n — 1 ) + e(n) 


• Described by: O = [1 1 1 2] => p = 15 


• Candidate terms: 

z(n) = 0o + 0iu{n) J r92u(n — l)+9^u 2 {n) + 6^u{n)u{n — l)-\-9^u 2 (n — 1) 
+ 9§z(n — l)+9^u{n)z{n — 1) + 9%u(n — 1 )z(n — 1) + 9$z 2 (n — 1) 
+ 9iQu(n)e{n — 1) + 9nu(n — 1 )e(n — 1) + 9i2z(n — 1 )e(n — 1) 

+ 9ize{n — l)+0i 4 e 2 (n — l)+e(n) 
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NARMAX Representation and Identification of Ankle Dynamics 
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Objectives 


• Theoretically derive a non-linear difference equation of a parallel pathway 
model of ankle dynamics 


• Show that the theoretical equation for this ankle model is of the NARMAX 
form 


• Investigate the suitability of NARMAX identification methods applied to 
ankle dynamics 
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Motivation 


• Better design of prosthetics 


• Characterisation of normal versus disease subjects 


• Early detection of neuromuscular disease 


• Design of robotic motor control systems 


• Astronaut health to mitigate muscle atrophy - develop optimal exercise 
regime 
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Parallel Pathway Model of Ankle Dynamics 


U{s 


Ankle 


Angle 


Is 2 + Bs + K 


Intrinsic Stiffness 
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Measured 


Torque^-^ Torque 

Y(s) 


Z(s 


Reflex Delay (Half-Wave Rectifier) Muscle Activation 


Reflex Stiffness 
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Parallel Pathway Model of Ankle Dynamics 
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Parallel Pathway Model of Ankle Dynamics 
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Parallel Pathway Model of Ankle Dynamics 
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Parallel Pathway Model of Ankle Dynamics 
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Parallel Pathway Model of Ankle Dynamics 
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Parallel Pathway Model of Ankle Dynamics 
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Theoretical Analysis: Linear Pathway 


U(s) 

Ankle Angle 


Is 2 + Bs + K 


Yl(s) 

► 

Linear Torque 


• Discrete-domain approximation to a derivative (Newton’s backwards for- 
mula) was used to approximate the intrinsic pathway dynamics 

d u(t ) u(n) — u(n — 1) 


where T = sampling rate and n = sampled data point index 
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Theoretical Analysis: Non-linear Pathway 



Reflex Delay (Half-Wave Rectifier) Muscle Activation 


• Continuous-time delay converted to discrete-time as r 
the continuous-time delay 


"A" 

T 


where A is 


• Half-wave rectifier approximated as a second-order static polynomial: cq + 

c\x(n) + C2X 2 (n) 


- Second-order fit accounted for over 98% of the output variance of static non-linearity 


Activation dynamics converted to discrete-time via the bilinear transform 

2 (z- 1 

T Vz + 1 


s — — 
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Theoretical Analysis: Overall Non-linear Model 

a fit . ~ ^ - , ,.5 *..^j -* 


• Collecting terms and combining yielded the overall non-linear model as 

z(n) = 60 + b\z{n — 1) + b 2 z(n — 2) + b^u(n) + b/\u(n — 1) + b$u(n — 2) 

+ bQu(n — 3) + bju{n — 4) + b%u(n — r) + bgu(n — r — 1) 

+ biou(n — r — 2) + b\\u{n — r — 3) + b\ 2 U 2 {n — r) + &i3it 2 (n — r — 
+ buu 2 (n — r — 2) + &i5?i 2 (ri — r — 3) + biQu(n — r)u(n — r — 1) 

+ — r — \)u(n — r — 2) + 6ig^(n — r — 2)u{n — r — 3) 

+ 6i 9 e(n — 1) + b 2 oe(n — 2) + e(n) 

• This is a NARMAX model since 

- (i) it includes input-output terms that are combinations of linear, non-linear integer 
powers and cross-products and 

- (ii) is linear-in-the-parameters 
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Validation of NARMAX Ankle Model 


• Simulate the NARMAX description model’s response and compare it to 
the continuous-time model’s response for typical parameter values found in 
experiments 


• Input sequence uniformly distributed, white, zero-mean, random input, ban- 
dlimited to 30 Hz 


• Operating range between ±0.40 rad 


• Sampling rate 200 Hz (T = 0.005 s) 
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Torque (Nm) Position (rad) 


Results 


Parallel Pathway and NARMAX Model Input 
0.4 

0.2 

0 


- 0.2 



Parallel Pathway and NARMAX Model Ouputs, 99.53%VAF 
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Identification of NARMAX Representation of Ankle Dynamics 

■ ~ ~ . ■ ■ ■ ------ ■ 


• Assess the utility of methods developed for identifying NARMAX models 


• Noise on the output enters the model as product terms with the system input 
and output 


• Ordinary least-squares yields biased estimate => Does not account for noise 


• One solution extended least-squares (ELS) =^> Bias addressed by modelling 
lagged errors 
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Identification of Full NARMAX Representation Yields Biased Estimates 


• Reduce number of terms 


• Not a general reduction of terms to describe the data but rather a minimi- 
sation of the number of regressors or degrees of freedom used to form the 
regressor matrix 


• Reduction provides a regressor matrix that is more stable in terms of invert- 
ibility since the coefficients will no longer be interrelated 


• Similar to reconditioning a matrix via normalisation 
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Reduced Model Representation 


• Recombining all terms according to coefficients of the static non-linearity 
yields an overall non-linear model represented by 12 terms as 

z(n) = 60 + b \ z(n — 1 ) + 62^(71 — 2 ) + 6377(77) + 6477(77 — 1 ) + 6577(77 — 2 ) 
+ bQu(n — 3) + 6777(77 — 4) + miv ( n ) + m2x{n) + bige(n — 1 ) 

+ 620^(77 — 2 ) + e(n) 

where 

77(77) = 77(77 — t ) + 77(77 — r — 1) — 77(77 — r — 2) — 77(77 — r — 3) 

and 

X(77) = 77 2 ( 77 — t ) + 3T7 2 (77 — T — 1) + 3t7 2 (t7 — T — 2) + 77 2 (t7 — T — 3) 

— 277(77 — t)t7(t7 — t — 1) — 477(77 — r — 1)77(77 — T — 2) 

— 277(77 — T — 2)77(77 — T — 3) 
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Experimental Data 


• Eight trials of experimental human ankle data analyzed from single subject 
with no history of neuromuscular disease 


• Pseudo-random binary sequences (PRBS) of 0.05 rad and 260 ms switching 
rate 


• Subject maintained a mean contraction of -5 Nm 
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Experimental Setup 
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Typical Position Input and Cross-Validated Torque Output 


PRBS Position Input 



Measured and Cross-validated Torque, 98.84%VAF 
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Cross-Validation % VAF for Each Trial 




Cross-Validation Fit 
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Conclusions 




• Non-linear difference equation describing parallel pathway model is a NAR- 
MAX model 


• Simulation results show that the NARMAX model matches the continuous- 
time response well 


• Identification methods developed for NARMAX models can be used to 
identify human ankle dynamics parametrically 


• Experimental data shows that estimated parameters explain the input-output 
data well 


• The NARMAX form is clearly amenable to the study of a wide range of 
biological systems 
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A Least-Squares Parameter Estimation Algorithm for Switched 
Hammerstein Systems With Applications to the VOR 
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Objectives 


• Can Vestibular Ocular Reflex (VOR) be described by a NARMAX model? 


• Develop identification technique to estimate parameters of non-linear hy- 
brid (switched) systems 


• Apply identification technique to experimental human VOR data 
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Motivation 


• Characterisation of normal versus disease subjects 


• Early detection and quantification of ocular disease and balance disorder 


• Design of robotic motor control systems 


• Pilot/astronaut health and safety 
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General Hammerstein Structure for Switched System 




Static NL M Linear Sub- System M 
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Extension of NARMAX Model for Non-Zero-Initial-State 


Zm(n) = f l m [u(ri ), . . . , u(n — n u )\ + £[z(n - 1), . . . , z(n - n z ), e{n - 1 ),..., e(n - n e )\ + 
zr m (n) = f l m [u(l ) : . . . , u(n — 1), u(n)] + £[£(n)] for m = 1, 2, . . . , M V finite M 
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Simulated VOR Data 


Head Velocity 



Eye Position 



0 1 2 3 4 5 6 

Time (s) 
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VOR Model 
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NARMAX Model of VOR 




• Dynamics converted to discrete-time via the bilinear transform 

• Collecting terms and combining yielded an overall non-linear model with 
two NARMAX sub-models: 


y{n) 
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f 2 / 1 H 

1 2/2 (n) 

Switch Position Si 
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+ Ph 
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i = !,-■■ ,q 
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Biased Parameter Estimate 




• Define: Z = + £ 

£ = Ai( 5 i(n — ti) + A 2A2 (n — <2) + • • • + X r S r (n — t r ) + e(n) 


£ 


e 


ELS 


yT E J Z ] = 0£ig + 




Bias: E (^ r ^)- 1 ^ r £| ^0 
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Modified Extended Least Squares 


• Modify ELS algorithm to correct bias due to initial conditions 


Include columns in the regressor matrix to account for the initial conditions 


$ = ht 


zue 




Extended parameter set based on this model formulation defined as: 
Qmf.t.r = (® T &) where Omf.t.r = \0 p 0$ 
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Algorithm 


• Segment all data according to mode 


• Concatenate data from each mode 


• Form regressor matrix for each mode 


• Add a column modelled as an impulse when new segment is used to form 
the regressor matrix 


• Estimate 6m els using minimisation 
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Experimental Data 


• Data analyzed from single patient with history of peripheral vestibular dis- 
ease - no function in the right inner ear following surgery 

- Often associated with large non-linearity in the VOR response and abnormally small 
slow-phase time constant 


• Experimental protocol used a sinusoidal rotation at 1/6 Hz, with a peak 
head velocity ~ 200 deg/s 

- Test 52 s, last 32 s recorded to measure VOR properties during sensory steady state 


• Signals sampled at 500 Hz then digitally low-pass filtered to 15 Hz to reduce 
high frequency content 
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Data Analysis 



Ui(s) = -^—-U(s) = 


S -U(s) U 2 (s) = ~^—U(s) = 


■U(s 


T\S + 1 s H- pi r 2 s+ 1 v/ S+P 2 

• Canal dynamics pre-process head velocity which serve as sensor for head movement 

• Sensory vestibular process well described by first-order, high-pass system which transmits 
signals to the central circuits acting as a switched system 

• Estimate T\^ associated with vestibular dynamics, assuming ri ;2 =1 -15 s in Is increments 

— Combination of vestibular time constant and switched system time constant which yielded highest cross- 
validation %QF was deemed the best-fit model 


• Compared the results of MELS to ELS (not modelling initial conditions) 
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Experimental Setup 


IYTiTi- ■ ■ 
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Experimental VOR Data 


Head Velocity 



Eye Position 
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Results: Cross- Validated Eye Position 


ELS Cross-Validated Slow-Phase Output: QF= 27.2% ELS Cross-Validated Fast-Phase Output: QF= 16.7% 
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MELS Cross-Validated Slow-Phase Output: QF= 98.5% 
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MELS Cross-Validated Fast-Phase Output: QF= 86.3% 
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Estimated Vestibular & Switched System Parameters 

'■ ! v-.i — ^ irr - , i ■ ia ■ ■ j. r .- ^ tM r, ifc E 




(a) Xq \ — 14.0 s, ti = 0.849 s 


(b) Tq 2 — 1-00 s, t 2 — 0.0366 s 



U(s) 


Xp,2S 
Xo, 2^+1 
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Dynamics 2 
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Third-Order MELS Fit to Fast-phase Nonlinearity 



MELS Static Nonlinearity Estimate 


*2 


J+3 
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MELS Fast-Phase 
Dynamics 


(c) Xqa — 10.0 s, ti = 0.915 s 


(d) to ,2 = 2.00 s, t 2 = 0.0831 s 
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■ jLJ&Xi 


Conclusions 


• MELS method provides accurate estimates of parameters since it takes ad- 
vantage of an entire data record even though the individual segments are 
short 


• Results may have a clinical significance in the analysis of ocular nystagmus 
of all types (pursuit, optokinetic, etc.) 


• Technique here allows greater insight into the functionality of various ocu- 
lar reflexes, by providing quantitative measures of both saccadic and slow 
ocular dynamics from a single experimental record 


• MELS method may be useful to estimate the coefficients of complex Ham- 
merstein structure switched systems in biology 
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A Suboptimal Bootstrap Method for Structure Detection of 
Non-linear Output-Error Models with Application to Human 

Ankle Dynamics 
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Objectives 


• Develop a structure detection method to determine a parsimonious model 
description which best fits the observed output 


• Gain insight into the underlying process describing ankle model 


• Assess whether morphological ankle model is accurate 
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Two Fundamental Approaches to Structure Detection 

’■ ir i.i.'mm ; — r -• i ■ BEi- ■ j . 7- ■-r.-ir.-_- ai-ifc A •* 


• Exhaustive search 

- Every possible subset of the full model is considered 

- Requires a large number of computations and known to not converge to true underlying 
system 


• Parameter variance 

- Parameter variance estimates computed from model residuals 

- Often inaccurate when number of candidate terms large 
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Bootstrap 


• Numerical procedure for estimating parameter statistics 


• Mild conditions on sample errors 

- Errors independent and identically distributed (i.i.d.) 

- Zero-mean 
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Hypothesis 


IYIYn ■ ■ 


The bootstrap might be useful for structure detection of 
over-parameterised non-linear models 


ICNPAA 2010 World Congress 


30 June - 3 July 2010 


IVTiTi- • 


Sao Jose dos Campos, Brazil 


Theoretical Analysis of Bootstrap 


• Bickel and Freedman (1981) analysed linear regression model where N and 
p both large 


• Showed for full p- dimensional distribution of least-squares estimates, boot- 
strap distribution will converge to true unknown distribution when 


7 = — ^ 0 « 0.1 


• Initially, p cannot change; accuracy of bootstrap estimate determined by 
data length, N, available for estimation 
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Data Length 




• Model order: O = [ 4 4 4 2] has p = 105 candidate terms 


• Data points: N = = 110, 250 


• Desirable to reduce number of candidate terms 
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Noise Model 


^hdi 


y(n) = 0\y(n — 4) + 92y 2 {n — 4) + 9^u 2 (n — 4 )y(n — 4) + 9^ u(n — 4) 


z(n) = y(n) + e{n 


l2 


z(ri) = 9i[z{n — 4) — e(n — 4)] + 92[z(n — 4) — e(n — 4) 

+ 9‘2 > u 2 {n — 4 )[z(n — 4) — e(n — 4)] + 9^u{n — 4) + e(n 


• Simple output additive noise can result in complex noise model 

• Number of noise terms often can be larger than process terms 

• Useful to avoid estimating noise model whilst yielding a unbiased estimate 

• Estimator? 
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Instrumental Variables 


• Based on selecting instrument matrix V which satisfies conditions 


lim — V T \I/ 

N—>oo N 


zu 


lim 

iV— >oo 




R; where R is nonsingular 
0 . 


- (i) the matrix product, R = \ T ^ ZU , has full rank, and 

- (ii) the errors have zero-mean and be uncorrelated with V 

• This ensures the estimate 

e IV = {v T y zu y l v T z 

is unbiased since the instrument matrix is not correlated with the errors 
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NARMAX with Linear Map 


IV gives unbiased estimates when noise represented as linear map 

z(n) = f l [ u (n), ■■■ , u(n — n u ), z(n - 1), • • • ,z(n-n 
+ C[e(n —!),••• , e(n — n e )} + e(n 


Assuming output additive noise 

z(n) = f l [u(n), ■■■ , u(n — n u )] + C[z(n - 1), • • • ,z(n- n z ), 
e(n —!),••• , e(n — n e )\ + e{n] 


Restricted class of NARMAX models 

- Blocked structured N-L models (static non-linearity followed by a causal, linear, time- 
invariant, dynamic system) such as Hammerstein models, bilinear models, etc. 
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Model Order & Candidate Terms 


• Redefine model order as O 


n u n z l } 


• Maximum number of candidate terms: 


i 

P = n z + ^pk + 1 
k = 1 

Pk-i{n u + k) 

Pk = 7 , Po = 1 


• Consider again earlier example: model of order: 0 = [4 44 2]=>p = 105 
candidate terms 

• Order redefined as O = [4 4 2] =^> p = 25 terms need be considered; a 
significant reduction 
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Parallel Pathway Model of Ankle Dynamics Revisited 


U{s 


Ankle 


Angle 


Is 2 + Bs + K 


Intrinsic Stiffness 


Ankle 


Velocity 
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Yl(s) 
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C Net 
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gu z 
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s 2 +2Ccjs+u; 2 



Measured 


Torque^-^ Torque 

Y(s) 


Z(s 


Reflex Delay (Half-Wave Rectifier) Muscle Activation 


Reflex Stiffness 
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Overall Non-linear Model of Ankle Dynamics 


z(n) = 60 + b\z(n — 1) + b 2 z(n — 2) + b^u{n) + b/±u(n — 1) + b$u(n — 2) 

+ bQu(n — 3) + b^u(n — 4) + b%u(n — r) + bgu(n — r — 1) 

+ biou(n — r — 2) + bnu(n — r — 3) + bi 2 U 2 (n — r) + b\^u 2 (n — r — 1) 
+ buu 2 (n — r — 2) + bi^u 2 (n — r — 3) + b\Qu(n — r)u(n — r — 1) 

+ bnu(n — t — 1 )u(n — r — 2) + b\%u{n — r — 2 )u(n — r — 3) 

+ bi 9 e(n - 1) + b 2 oe(n - 2) + e(n) 

• Reflex delay 50-100 ms; corresponding to a discrete-time delay r = 

- Maximum model order O = [8 2 2 2] 

• Gives full model description with 105 candidate terms 

- True system described by 21 parameters 

• Computationally intractable combinatorial optimisation problem 
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Overall Non-linear Model of Ankle Dynamics 


z(n) = 60 + b\z(n — 1 ) + b^zin — 2 ) + b^u{n) + b/^uin — 1 ) + b$u(n — 2 ) 

+ bQu(n — 3) + b^u{n — 4) + b%u(n — r) + bgu(n — r — 1) 

+ biou(n — t — 2) + b\\u{n — r — 3) + bi 2 U 2 (n — r) + bisu 2 (n — r — 1) 

+ buu 2 (n — t — 2) + bi^u 2 (n — r — 3) + biQu(n — r)u(n — r — 1) 

+ bi 7 u(n — t — 1 )u(n — r — 2) + bisu(n — t — 2 )u(n — r — 3) 

+ b ld e(n — 1 ) + 6 20 e(n — 2 -)-+ e(ri) 

• Reflex delay 50-100 ms; corresponding to a discrete-time delay r = 

- Maximum model order O — [ 8 2 2 2] 

• Use a priori knowledge to eliminate unrealistic candidate terms, increasing 
computational efficiency 

- Eliminate all nonlinear and cross-terms associated with discrete-time delays between 

r = [0 2 ] 

• Gives full model description with 33 candidate terms 

- True system described by 19 parameters 

• Computationally tractable combinatorial optimisation problem 
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Assumptions 


• System order assumed known 

- Small signal perturbations to estimate a LTI model order 


• SOBSD algorithm limits possible class of models but it is reasonable to 
consider this limited set since the physical model suggests that the true 
system lies within this class 


• Structure detection provide useful process insights that can be used in sub- 
sequent development or refinement of physical models 
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Experimental Data 


• Data analysed two data sets from a subject with no history of neuromuscular 
disease 

• PRBS input with 0.06 rad peak-to-peak amplitude and switching rate of 125 
ms and 260 ms whilst subject maintained a mean contraction of -5 Nm 

• Input-output recorded for 50 seconds 

• Measured data sampled at 500 Hz, for estimation data decimated to 50 Hz 

• Estimation data N e = 2,000, validation data N v = 1,000 and B = 100 
bootstrap replications to estimate parameter statistics 
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Torque (Nm) 


Results of Identifying Experimental Human Ankle Data Set 1 


Position Input 



Measured and Cross-validated Torque, 95.47%QF 



Time (s) 
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Results of Identifying Experimental Human Ankle Data Set 2 


Position Input 



Measured and Cross-validated Torque, 96.01 %QF 
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Conclusions 


• SOBSD algorithm provides a novel approach to the model selection prob- 
lem without computing a noise model and, hence, reducing the dimension- 
ality of this ill-posed combinatorial optimisation problem 

• Reduction in dimensionality comes at cost of limited model structures that 
can be considered 

• Study illustrates usefulness of structure detection as an approach to validate 
morphological models via analysis of input-output data 

• Results show identified model structure matches the theoretically expected 
structure well 

• Indicates morphological modeling studies may be accurate for this model 
describing ankle dynamics 
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